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ROSENBLATT DISTRIBUTION SUBORDINATED TO GAUSSIAN 
RANDOM FIELDS WITH LONG-RANGE DEPENDENCE 

N.N. LEONENKO, M.D. RUIZ-MEDINA, AND M.S. TAQQU 


Abstract. The Karhunen-Loeve expansion and the Fredholm determinant formula are used, 
to derive an asymptotic Rosenblatt-type distribution of a sequence of integrals of quadratic 
functions of Gaussian stationary random fields on R d displaying long-range dependence. This 
distribution reduces to the usual Rosenblatt distribution when dm 1. Several properties of this 
new distribution are obtained. Specifically, its series representation, in terms of independent chi- 
squared random variables, is established. Its Levy-Khintchine representation, and membership 
to the Thorin subclass of self-decomposable distributions are obtained as well. The existence 
and boundedness of its probability density then follow as a direct consequence. 
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1. Introduction 

The aim of this paper is to derive and study the properties of the limit distribution, as 
T —> oo, of the random integral 

(1.1) S T = 2- [ (U 2 (x) - 1 )dx, 

«r Jd(t) 

where the normalizing function dx is given by 

(1.2) d T = T d ~ a C(T), 0 < a < d/2, 
with C being a positive slowly varying function at infinity, that is 

(1.3) lim £(T||x||)/£(T) = 1, 

T-> oo 

for every ||x|| > 0, and D(T ) C M. d denotes a homothetic transformation of a set D C M d , with 
center at the point 0 € D, and coefficient or scale factor T > 0. In the subsequent development, 
D is assumed to be a regular bounded domain, whose interior has positive Lebesgue measure, 
and with boundary having null Lebesgue measure. Here, {T(x), x G is a zero-mean Gauss¬ 
ian homogeneous and isotropic random field with values in R, displaying long-range dependence. 
That is, Y is assumed to satisfy the following condition: 

Condition Al. The random held {Y (x), x G M d } is a measurable zero-mean Gaussian homo¬ 
geneous and isotropic mean-square continuous random held on a probability space (£l,A,P), 
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with E7 2 (x) = 1, for all x £ M d , and correlation function E[F(x)K(y)] = R(||x — y||) of the 
form: 

(1.4) B( || Z ||) = ^M), z€M rf , 0 < a < d/2. 

\\ z \\ 

From Condition Al, the correlation B of Y is a continuous function of r = ||z||. It then follows 
that C{r) = 0(r a ), r —> 0. Note that the covariance function 

(1-5) B(\\z\\) = (1 + || Z ||£) 7 ’ 0< ^ 2 > 7>0, 

is a particular case of the family of covariance functions (11.41) studied here with a = /Ty, and 

(1.6) £(M) = l|z||' 57 /(i + INI 5 ) 1 . 

The limit random variable of (11.11) will be denoted as Soo- The distribution of will be 
referred to as the Rosenblatt-type distribution, or sometimes simply as the Rosenblatt distribu¬ 
tion because this is how it is known in the case d = 1. In that case, a discretized version in time 
of the integral (11.11) first appears in [35], and the limit functional version is considered in [32] 
in the form of the Rosenblatt process. In this classical setting, the limit of is represented 
by a double Wiener-ltd stochastic integral (see 02; 03]). Other relevant references include, for 
example, im [is], in], m, [36] . to mention just a few. The general approach considered 
here for deriving the weak-convergence to the Rosenblatt distribution is inspired by [42] . which 
is based on the convergence of characteristic functions. This approach has also been used, re¬ 
cently, in [23], to study the characteristic functions of quadratic forms of strongly-correlated 
Gaussian random variables sequences. 

We suppose here d > 2, and thus consider integrals of quadratic functions of long-range 
dependence zero-mean Gaussian stationary random fields. We pursue, however, a different 
methodology than in the case d = 1, which was based on the discretization of the parameter 
space. A direct extension of these techniques is not available when d > 2. Instead of discretizing 
the parameter space of the random field, we focus on the characteristic function for quadratic 
forms for Hilbert-valued Gaussian random variables (see, for example, m, and take advantage 
of functional analytical tools, like the Karhunen-Loeve expansion and the Fredholm determi¬ 
nant formula, to obtain the convergence in distribution to a limit random variable S 0 0 with 
Rosenblatt-type distribution. 

The double Wiener-Ito stochastic integral representation of Soo in the spectral domain leads 
to its series expansion in terms of independent chi-squared random variables, weighted by the 
eigenvalues of the integral operator introduced in equation (13.11) below. The asymptotics of these 
eigenvalues is given in Corollary [2j The infinitely divisible property of Soo is then obtained as 
a direct consequence of the previous results derived, in relation to the series expansion of Soo, 
and the asymptotic properties of the eigenvalues. We also prove that the distribution of Soo is 
self-decomposable, and that it belongs, in particular, to the Thorin subclass. The existence and 
boundedness of the probability density of Soo then follows. 

The outline of the paper is now described. In Section [2] we recall the Karhunen-Loeve ex¬ 
pansion, introduce the Fredholm determinant formula, and use the referred tools to obtain the 
characteristic function of (11.11) . In Section [3] we prove the weak convergence of (11.11) to the ran¬ 
dom variable Soo with a Rosenblatt-type distribution. The double Wiener-Ito stochastic integral 
representation of Soo, its series expansion in terms of independent chi-square random variables, 
and the asymptotics of the involved eigenvalues are established in Section 0J These results are 
applied in Section [5] to derive some properties of the Rosenblatt distribution, e.g., infinitely 
divisible property, self-decomposability, and, in particular, the membership to the Thorin sub¬ 
class. Appendices A-C provide some auxiliary results and the proofs of some propositions and 
corollaries. 
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In this paper we consider the case of real-valued random fields. In what follows we use the 
symbols C. Cq, M\. M 2 , etc., to denote constants. The same symbol may be used for different 
constants appearing in the text. 


2. KARHUNEN-LOEVE EXPANSION AND RELATED RESULTS 

This section introduces some preliminary definitions, assumptions and lemmas hereafter used 
in the derivation of the main results of this paper. We start with the Karhunen-Loeve Theorem 
for a zero-mean second-order random field {Y(x), x € K C M^}, with continuous covariance 
function £?o(x,y) = E[Y(x)Y(y)], (x,y) E K x K C x M. d , defined on a compact set K of 
(see Section 3.2 in Cl)- This theorem provides the following orthogonal expansion of the 
random field Y : 

OO 

y ( x ) = Y x e K - 

3 = 1 

Afc^x) = / B 0 {x,y)(j) k (y)dy, k E N*, (fa,fa) L 2 (K) = b J € N*, 

J K 

(2.1) 

where r] k = ~^= j K Y (x)0 fc (x)dx, for each & > 1, and the convergence holds in the L 2 (fi, A, P ) 
sense. The eigenvalues of Bq are considered to be arranged in decreasing order of magnitude, 
that is, Ai > A 2 > • • • > Afc_i > A&>.... The orthonormality of the eigenfunctions cj)j , j € N*, 
leads to the uncorrelation of the random variables rjj, j E N*, with variance one, since 

E i r 1jVk]= / - B o(x,y)0 j (y)^ fc (x)dydx = Xj / <^-(x)<^(x)dx = Aj<5j jfc , 

JKJK JK 

with 5 denoting the Kronecker delta function. In the Gaussian case, they are independent. 

For each T > 0, let us fix some notation related to the Karhunen-Loeve expansion of the re¬ 
striction to the set D(T ) of Gaussian random field Y. with covariance function HI. 41) . By Ry,d(t ) 
we denote the covariance operator of Y with covariance kernel -Bo,t(x, y) = E[Y(x)Y(y)], 
x, y € D(T), which, as an operator from L 2 (D(T )) onto L 2 (D(T)), satisfies 

Ry,d(T)(4>i,t)(a) = / -Bo,T(x,y)0 /T (y)dy = Ai,r(fiy,D(r))('i,r( x )> 1 £ N*, 

Jd{t) 

where, in the following, by Xk(A) we will denote the /cth eigenvalue of the operator A. In partic¬ 
ular, {Afc i 7’(i2y,D(T))}fcLi and {^kT\V= 1 respectively denote the eigenvalues and eigenfunctions 
of Ry,d(t) , for each T > 0. Note that, as commented, i?o,T refers to the covariance function of 
{Y(x), x € D(T)} as a function of (x, y) E D(T) x D(T), which, under Condition Al, defines 
a non-negative, symmetric and continuous kernel on D(T), satisfying the conditions assumed 
in Mercer’s Theorem. Hence, the Karhunen-Loeve expansion of random field Y holds on D(T), 
and its covariance kernel L?o,t also admits the series representation 

OO 

(2-2) B 0}T (x,y) = Y Aj,T( J Ry,D(T))^,r( x )^,T(y) ) x -ye D(T), 

3 = 1 

where the convergence is absolute and uniform (see, for example, [1], pp.70-74). The orthonor¬ 
mality of the eigenfunctions {fa t\iLi yields 


1 

dx 



K 2 (x)dx = — Y X 3,t( r y,d(t))vIt- 
1 i=i 


(2.3) 
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In the derivation of the limit characteristic function of dm we will use the Fredholm de¬ 
terminant formula of a trace operator. Recall first that a positive operator A on a separable 
Hilbert space H is a trace operator if 

(2.4) MU, = TV(.4) ee £ {(A'A^ t , Vk ) < co, 

k 


where A* denotes the adjoint of A and is an orthonormal basis of the Hilbert space H (see 
[53] . pp. 207-209). A sufficient condition for a compact and self-adjoint operator A to belong 
to the trace class is J2T=i ^k(A) < oo. For each finite T > 0, the operator R-y.d(t) is in the 
trace class, since from equation (12.21) . applying the orthonormality of the eigenfunction system 
3 € N*}, and keeping in mind that Hq.t(O) = 1, we have 


OO r. 

(2.5) Tr(R Y ,D(T)) = Y ^3 ,t(Ry,d(t)) = / 

j= i Jd(t) 


R 0 ,T(x,x)dx = 


Id(t) 


dx = T d \D\ < oo, 


where \D\ denotes the Lebesgue measure of the compact set D. Note that the class of compact 
and self-adjoint operators contains the class of trace and self-adjoint operators. Hence, under 
Condition Al, from equation (12.51) . the restriction of Y to D{T ) admits a Karhunen-Loeve 
expansion, convergent in the mean-square sense (i.e., in the L 2 (Q,A.,P)- sense), for any T > 0, 
and for an arbitrary regular bounded domain D. Furthermore, for any k > 1, 

(2.6) RY, D (T)fW= f <?(x,y)/(y)dy, / g l 2 (d(t)), 

Jd(t) 

where B^} denotes 

B o!t ( x > y ) = b o,t ( x , y ), k = 1 , 


(2.7) B*^\x,y) = [ b * ( t 1) (x,z)R 0)T (z,y)dz, k = 2,3,.... 

Jd(t) 

From equations ()2.2|) and (12.71) . applying the orthonormality of 4>j,Ti 3 S N*, one can obtain 

OO « 

(2.8) Tr(Ry D(T) ) = VA j T ( R Y,D(T)) = / Ho ( r(x,x)dx < oo, k € N*, 

“7 Jem 


i=i 


since, for every A: > 1, |(-R^z)(T’))I A M\\y.{R. YD ( T ))\ k = M| X k {Ry d(T))I, f° r some positive 
constant M. In particular, in the homogeneous random field case, 


Br ( B Y,D(T)) 


(2.9) 


Y ^J,t( B Y,D(T )) 

3 = 1 



B, 


*(k) 
0 ,T 


(x k ,x k )dx k 



and, in the homogeneous and isotropic case, for k = 2, 


(2.10) Tr(Ry D(T )) — y; \\ t (Ry,d(t)) 

3 =1 


f r £ 2 (||x-y|| 
Id{t)Jd{t) ll x — y|| 2a 


-dydx. 


The following definition introduces the Fredholm determinant of an operator A, as a complex¬ 
valued function which generalizes the determinant of a matrix. 
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Definition 1 . (see, for example, [39] . Chapter 5, pp.47-48, equation (5.12)) Let A be a trace 
operator on a separable Hilbert space H. The Fredholm determinant of A is 


( 2 . 11 ) 


/ TSL TrA k 

T>(u) = det(I — uiA) = exp | — ^2 —-— 


A I = 


opI-EElwii 


k =1 


k= 1 1=1 


for u € C, and |cu|||Al||i < 1. Note that ||^4 m ||i < for A being a trace operator. 

Lemma 1. Let {Y(x), x E D C M d } be an integrable and continuous, in the mean-square 
sense, zero-mean, Gaussian random field, on a bounded regular domain D Cff* containing the 
point zero. Then, the following identity holds: 


E 


exp 


i£ [ Y 2 (x)d> 

Jd 


( 2 . 12 ) 


[](1 - 2X j (R Y , D )iO~ 1/2 = 2iO)~ 1/2 


for ||i?y,u ||i 12i^| < 1, as given in Definition Ql 

Proof. The covariance operator Ry,d of Y, acting on the space L 2 {D), is in the trace class. From 
Definition [H the following identities hold: 


E 

exp(i£ f Y 2 (x)dx) 

= E 

exp ( HZWRrJfg ] 


\ Jd J. 


L v U )\ 


I] E [exp (ifXjiRY^ri 2 )] = [](1 - 2X j {R Y , D )i(f)- 1/2 = (Vm))~ 1/2 

3 =1 3=1 


exp 


m= 1 


m 


“P ( \ E ^Tr(iJ^ D; , , 


2 *—' m 

m=l 


(2.13) 

where the last two identities in equation (12.131) are finite for |£| < from the Fredholm 

determinant formula (12.111) . Note that 

oo oo 

(2.14) Tr(R YT) ) = ^2 < X'f x ( R y ,d ) ^2 ^ 3 (Ry,d) = A™ 1 (Ry,d)\\R y ,v\\i < oo. 

3=1 3=1 

□ 


Remark 1. Similarly to equation H2.12\) . one can obtain the following identities, which will be 
used in the subsequent development: For a homothetic transformation D ( T) of D C with 
center at the point 0 € D, and coefficient T > 0, 

OO 

= l\(l - 2X j , T (R YM T)m~ 1/2 = <Pt{ 2iO)" 1/2 

\ m= 1 / 

where X i,t(R Y ,d(t)) > M,t(Ry,d{t)) >■■■> Aj ,t{Ry,d(t )) > • • •, with, as before, {X j, T (R Y ,D(T)), 3 € 
N*} denoting the system of eigenvalues of the covariance operator R y ,d(t) °fY, as an operator 
from L 2 (D(T )) onto L 2 (D{T)). The last identity in equation 112.15 )) holds for H-Ry^T) 111|2i£| < 

1, i.e., for Tr(R Y> D(T))\2i£\ = T d \D\\2i£\ < 1, or equivalently for |f| < 2T ] |D| ■ 
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3. Weak convergence of the random integral St 

This section provides the weak convergence of the random integral m to a Rosenblatt-type 
distribution, in Theorem [2j This results is based on the asymptotic behavior of the eigenvalues 
of the integral operator IC a (see Theorem |T] below) 


(3.1) 


£<*(/) ( x ) = 


Id ll x -y|l 


;f(y)dy, v/ € Supp (JC a ), 0 <a <d, 


with Supp(A) denoting the support of operator A. Operator (13.11) can be related with the Riesz 
potential (—A)~d/ 2 of order /3, 0 < (3 < d, on R d , formally defined as (see [H] . p.117) 

(3.2) (-A)-W 2 (/)( x ) = J-/ ||x — y\\~ d+ ^ f(y)dy, 

T\P) Jn. d 

where (—A) denotes the negative Laplacian operator, and 

Ti d / 2 2 l3 r((3/2) 1 


(3.3) 


l(P) = 


d—/3 

2 


c(d,/3)’ 


0 < (3 < d. 


Indeed, except a constant, the function (l/||x — y|| a ) in equation (13.11) defines the kernel of the 
Riesz potential (—A)(" _d )/ 2 of order [3 = (d — a), for 0 < a < d. Similarly, (l/||x — y|| 2a ) is 
the kernel of the Riesz potential (—A)“ _rf / 2 of order (3 = (d — 2a) on R d , for 0 < a < d/2. 

Recall that the Schwartz space S(R d ) is the space of of infinitely differentiable functions on 
M d , whose derivatives remain bounded when multiplied by polynomials, i.e., whose derivatives 
are rapidly decreasing. Particularly, Cq°(D ) C 5(M d ), with Cq°(D) denoting the infinitely 
differentiable functions with compact support contained in D. 

The Fourier transform of the Riesz potential is understood in the weak sense, considering the 
space S(R d ). The following lemma provides such a transform (see Lemma 1 of [41), p. 1 IT): 


Lemma 2. Let us consider 0 < (3 < d. 

(i) The Fourier transform of the function ||z|| — is 7(/3)||z|| - ^, in the sense that 
(3.4) f \\z\\~ d+l 3 if(z)dz = f 7(/3)||z|| _/3 J r (?/;)(z)(lz, Vi/’ € <S( R d ), 

J R d JR d 


where 


F(if)(z) = / exp (— i (x, z)) ip(x)dx 

J R d 


denotes the Fourier transform ofi/>. 

(ii) The identity T ((—A)~^ 2 (/)) (z) = ||z|| _|5 J r (/)( Z ) holds in the sense that 

(3.5) f (-A)' /3/2 (/)(x) 5 (x)dx = — irj f J 7 (/)(x)||x||^J r (c/)(x)dx, V/,c/€5(M d ). 

JR d (i7T) a J R d 

In particular, the following convolution formula is obtained by iteration of (13.51) using (13.21) : 


!m. d \ j(P) Jm 

= [ (-A) 

JR d 


x - y 


i -d+p 


1 


7(/3) J 

(-A)-^ 2 (f)] (x) ^T)d: 


y - z 


| -d+0 


f(z)dz 


dy 3(x)dx 


(2tt) 


j-((_A)-/3/ 2 (/))(x) 


X 


x 


-^( 5 )(x)d: 


x 
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J-(/)(x)||x||^||x||-^)(x)d: 


(2vr) d JK „ 

Kz f jr (/)( x )ii x ir 2/3 -Hff)(x)d> 

" ) J R d 


X 


(2vr 

= [ (-^) _/3 (/)( x )ff( x )^ x > V/,5 £ S(R d ), 0 < /3 < d/2, 

(3.6) 

where we have used that if / € 5(R rf ), then (— A)~P/ 2 (f) € 5(R d ). From equation (|3.6I) . and 
Lemma [2](i) , 


1 


/Rrf 7( 2 /3) 


|-^/(z)dz= / ||z||- 2/3 J-(/)(z)(/z 


1 


/Rd [7(/3)P 


l z - y| 


" d+/3 iiyir d+/3 d> 


/(z)dz, V/ e S(R d ), 0 < p < d/2. 


(3.7) 


Let us now consider on the space of infinitely differentiable functions with compact support 
contained in D, Cfi°(D) C 5(R d ), the norm 


II/IIU— = = ((-A r-^(f),f) LHD) 

= / R ,(-™^ = 1 1 

(3.8) =7^Td[ \Hf)W\ 2 \\M\- {d - 2a) dX, V/ € C™(D), 0 < a < d/2. 

(2vr) a 7 Rd 

The associated inner product is given by 

(3 ' 9) = / D 7^/ D 1^7F /(y) ^ y<ix ' 

for all f,g € Cq°(D). The closure of Cq°(D) with the norm || ■ ||(_A)“- d / 2 > introduced in (13.81) . 

defines a Hilbert space, which will be denoted as H. 2 a-d = Cq°(-D/ K-^) a ~ d/2 _ 

Remark 2. For a bounded open domain D, from Proposition 2.2. in [9], with D = n — 1, 
p = q = 2, and s = 0 (hence, A pq (D) = A^iD) = L 2 (D), where, as usual, L 2 (D) denotes the 
space of square integrable functions on D), we have 

(3.10) C^%D) IMIi2(Kd) = L 2 (D), 

(see also )45j . for the case of regular bounded open domains with C°°— boundaries). In addition, 
for all f € Cq°(D), by definition of the norm iS.8\) . 

Il/Il(-A)“—d/2 < C'||/||l / 2( R d), 

that is, all convergent sequences of Cq°(D ) in the L 2 ( R rf ) norm are also convergent in the 
hi. 2 a—d norm. Hence, the closure ofCo°(D), with respect to the norm || • || L 2 ( R d), is included in 
the closure ofC(/°(D), with respect to the norm || ■ ||(_A)“- d / 2 - Therefore, from equation \3.10 1) . 

l 2 ( d ) = c ^( d ) Ul2 ^ c qf ( L >) IMI (- A >“- d/2 = n 2 a - d - 


(3.11) 
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The asymptotic order of the eigenvalues of operator JC a on L 2 (D), in the case d > 2, are 
given in the next result (see, for example, [47], [49] and [SQ, p.197). (See also [13] and [48], for 
the case d = 1). 


Theorem 1. Let us consider the integral operator IC a introduced in equation Id. II) as an operator 
on the space L 2 (D ), with D being a bounded open domain of R rf . The following asymptotic is 
satisfied by the eigenvalues A/ c (/C Q ,), k > 1, of operator KL a : 


(3.12) 


r A fc(/C Q ) 

iim -—- 7-j—r 

fc—Aoo k-( d ~ a )/ d 


= T(d,a)\D\( d - a V d , 


where \D\ denotes, as before, the Lebesgue measure of domain D , and 

(d—a)/d p ^ d—a ^ 


(3.13) 


c{d, a) = tt“ /2 


T(f) [r (|)] 


( d—ot)/d ' 


Proof. We apply the results derived in [39], on the asymptotic behavior of the eigenvalues 
associated with certain class of integral equations. Specifically, the following integral equation 
is considered in that paper: 


(3.14) 


J I/ 1/2 (x)/c(x - y)V 1/2 (y)f(y)dy = A/(x), 


where k is an integrable function over a Euclidean space Ed of dimension d, having positive 
Fourier transform, and where V is a bounded non-negative function with bounded support. In 
particular, [49j considers the case where Ed = M d , V is the indicator function of a bounded 
domain D C R d , and k {||x — y||) = ||x — y|| a , for a > —d, and a / 0,2,4,.... Function k 
coincides in M d \ D with a function whose Fourier transform /(£) is asymptotically equal to 

r ( d ~ a 1 

r\d—ot d/2i \ 2 / | £ |— d-\-ot 

r (§) 4 

(see also the right-hand side of equation (13.41) for /3 = d — a, with 0 < a < d). For a > 
—d, a ^ 0, 2,4,... , the following asymptotic of the eigenvalues of the integral operator with 
kernel k{ ||x — y||) = ||x — y||“ is given in equation (2) in [35] : 


d+a. 

(3.15) Afc ~ 7 T~r 2 ( - 


r(fe) 


dj r(^) [r(f)] <d+ “ )/d 


[V(x)] i,id+a> dx 


(d+a)/d 


d+a)/d 


with 




d/ (d—a) , _ 


dx = \D\. 


Since function k in [39] coincides with the kernel of the integral operator IC a in equation (13.11) . 
for a € (—d, 0), equation ([3.15D then leads to the following asymptotic of the eigenvalues of JC a : 

r r . -| (d-a)/d 

A k (Ka) ~ »“ /2 ^ 


r(^) 

d > r(f)[r(|)]" i -“ )/ ' i 


I [V(x)| 

Jwi d 


d/(d—a) 


dx 


fo-{d-ot)/d' 


□ 


Remark 3. Similar results to those ones presented in Theorem 3.2 of [48] can be derived for the 
spectral zeta function of the Dirichlet Laplacian on a bounded closed multidimensional interval 
ofR d (see also m, for the case of d = 1). For a continuous function of the negative Dirichlet 
Laplacian, the explicit computation of its trace cannot always be obtained in a general regular 
bounded open domain of W l . Specifically, the knowledge of the eigenvalues is guaranteed for 
highly symmetric regions like the the sphere, or regions bounded by parallel planes (see, for 
example, [so]; m; mv- In particular, for the torus T 2 in M 2 , the Spectral Zeta Function can 
be explicitly computed (see, for example, [6], Chapter 1, equation (1-49), pp. 28-29). 
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The following condition is assumed to be satisfied by the slowly varying function C, in Theorem 
[2] below. 


Condition A2. For every m > 2, there exists a constant C > 0 such that 
£(T||xi — x 2 1|) £(T||x 2 - x 3 1|) £(T||x m -xi 


(3.16) 


Id ^ Jd £( t )II x i - x 2 || a T(r)||x 2 - x 3 ||“ C{T) ||x m - xi| 

dxidx -2 ■ ■ ■ dx m 


dxidx2 • • • dx m < 


< C 


'D 


-M- / 

J L 


D || x l - x 2|rH x 2 - x 31 


x m - X 1 


Note that Condition A2 is satisfied by slowly varying functions such that 


(3.17) 


sup 

T,xi,X2G-D 


£(T|| X i - x- 2 ||) 
£(T) 


<C 0 , 


for 0 < Co < 1. This condition holds for bounded slowly varying functions as in (II. 61) . in the 
case where D C Si(0), with £q(0) = {x € M d , ||x|| < 1}. 

For the derivation of the limit distribution, when T — > 00 , of the functional (EH), we first 
compute its variance, in terms of H 2 , the Hermite polynomial of order 2. It is well-known that 
Hermite polynomials form a complete orthogonal system of the Hilbert space L 2 (M, ip(u)du), 
the space of square integrable functions with respect to the standard normal density ip. They 
are defined as follows: 

Hk(u) = (-l) fe e^J^e _ ^, A: = 0,1,.... 

In particular, for a zero-mean Gaussian random field Y, for k > 1, 


(3.18) E H k (Y(x)) = 0, E (H k (Y(x)) H m (Y{ y))) = 5 m , k ml (E[y(x)T(y)]) m 
(see, for example, [33|). 

We use some ideas from m , Sections 1.4, 1.5 and 2.1). Consider the uniform distribution 
on D(T) with the density: 

(3.19) Po( T )(x)=T d \D\ 1 I xe D('r); x € ! d , 

where ^ x ^o(T) denotes the indicator function of set D(T). 

Let U and V be two independent and uniformly distributed inside the set D(T ) random 
vectors. We denote V’_d(t)(/ 9 )> the density of the Euclidean distance ||U — V||. Note that 
V’d(t) (p) = 0, if p > diain (D(T )), and 4>D(i){p ) is bounded, where diam (D(T)) is the di¬ 
ameter of the set D(T). 

Using the above notation, we obtain 


/ / 

Jd(t) Jl 


D(T) JD(T) 


G(||x-y||)dxdy 


\D(T)\ 2 E [g (IIu — V||)] 


(3.20) 


/•diam(D(T)) 

D\ 2 T 2d / G(p)iP D , T) (p)dp, 


J 0 

for any Borel function G such that the Lebesgue integral (13.201) exists. In particular, under 
Conditions Al—A2 for 0 < a < d/2, and T —>• 00 , we obtain 


ct 2 (T) = Var 


lD(T) 


H 2 (Y (x))dx 


= 2 


Id{t)Jd(t) ii x — y 




(3.21) 


= 2\\D\ T' 


2 qn2d 


l 


diam(D(T)) 


^ (P)P a *pD(T)(p)dp- 
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In equation (|3.21l) . consider the change of variable u = p/T. Applying the consistency of the 
uniform distribution with a homothetic transformation, and the asymptotic properties of slowly 
varying functions (see Theorem 2.7 in (38]) we get 

/■diam(D) 

2rp2d—2a / U - 2a C 2 {uT)^ D (u)du 


a{T) = 2\D\ Z T 


(3.22) 

where, by (13.201) . 


_ \p\2rp2d-2ot r 2 


C z (T)[ad{D)Y (1 + o(l)), 0 < a < d/2, T —>• oo, 


(3.23) 


ad(D) = 


■■L 


diarn(D) 


1 1/2 


u 2a if D (u)du 



dx.dy 


1 1/2 


dJd ii x - y 


1 2a 


More details, including properties of slowly varying functions, can be found in [3]. 
If D is the ball Bt{ 0) = (x € : ||x|| < T}, then (see [17], Lemma 1.4.2) 


Yb t (P){p) ~ T d/ i -(^) 2 ( 


yyy ) dp, 0 < p < 2T, 


(3.24) 
where 

(3.25) 

r (p) r (?) Jo 

is the incomplete beta function. In this case, one can show (see Lemma 2.1.3 in H3) 

2d-2a+2 n d-l/2Y ^d = 2a± 1) 


4(P»g) = f fP 1 ( 1 ~ t ) q ldt , /i €[0,1], p>0, q > 0, 

r(p)r(g) j o 


(3.26) a d (fli(0)) = 

For d = 1, D = [0,1], 

oi([0,l]) = 2 


(d — 2 a)r (|) T(d — a + 1) 

1 rl dxdy 



1 


, 0 < a < 1/2. 


/o Jo \x~y\ 2a (1 - a)(l - 2a) : 

Let us now consider domain D = [0, l\] x • • • x [0, If C M d , li > 0, i = 1,..., d. The Dirich- 
let negative Laplacian operator on such a domain has eigenvectors {4>k}k>i and eigenvalues 
{A fc (-A)} fc >i given by 

d 


4( x ) = II sin 


7T kiXi 


i —1 

d 7 r 2 Z ’ 2 

•M-a) = 

j=l 




, X = (xi,... ,x d ) G [0, Zi] x • • • x [0, Z rf ] , ki > 1, i = l,...,d 


, k = (ki,...,k d ), ki> 1, i = l,...,d. 


(3.27) 

The norm ||/C^,||/:(^ 2 (z?)) °f in B(L 2 (D)), the space of bounded linear operators on L 2 (D), 
is given by the supremum of its eigenvalues. From Theorem 4.5(h) in [10], 

a—d 


(3.28) 


A fc (/C*) < 2 


£ 


7 r 


! k 2 


If 

_i=l 1 


, d — a € (0,1). 


Note that Theorem 4.5(h) in [TO] holds for a bounded convex domain in M. d . From equation 


\\^-a\\c{L 2 {D)) < 1- 


(I3.28P , if min{/i ,■■■,ld} < \/ §7T, 
(3.29) 


Theorem 2. Let D be a regular bounded open domain. Assume that Conditions A1 and A2 
are satisfied. The following assertions then hold: 
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(i) As T —y oo, the functional St in (EP converges in distribution sense to a zero-mean 
random variable Soo, with characteristic function given by 


(3.30) 'ip(z) = E [exp(izS' 00 )] = exp f i ^ c m J , 

\ m=2 171 ) 


z e 


where, for m > 2, 
c m = Tr(K%) 


L 


(3.31) 

(ii) The functional 


D M Jd II x i - x 2|r 11 x 2 - x 31 


x m - X 1 


-dx i... dx„,. 




£{T)T d - 


[ G{Y(x))dx-C^T d \D\ 

JD(T) 


converges in distribution sense, as T — y oo, to the random variable \C^ S^, with S^ 
having characteristic function \3.30\) . and with G £ L 2 (R, <p(x)dx) having Hermite rank 
m = 2. Here, 

Cq = [ G(u)Ho(u)ip(u)du = E[G(y(x))] 

Jr 

C 2 = [ G(u)H- 2 (u)(p(u)du, 

Jr 

respectively denote the 0th and 2th Hermite coefficients of the function G. 

Remark 4. Note that Condition A2 is satisfied by the slowly varying function m with 
(7 = 1, for D = £>i(0) = {x : ||x|| < 1}. 

Proof. We first prove (i). Since EY 2 (x) = 1, 


[ dx= I E [y 2 (x)] dx = £ 

[ y 2 (x)dx 

/ D(T) J D(T) 

J D(T) 


OO OO 

= X] ^hT( R Y,D(T)) Er i] = X] ^3,t( R Y,D(T))- 

3 = 1 j =1 


From Definition [H Lemma [TJ and Remark [L] one has 




E 


exp 





lz Yl’jLl ^ j,T{Ry,D(T )) 
dT 



1 - 2\z 


^j,T(Ry,D(T )) 
dx 


- 1/2 



S£ll ^ 3,t{ R Y,D(T )) 
3t 

iz J2°jL 1 ^j,T{ R Y,D(T )) 

dx 
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= exp 


-\?',T(-Ry,D(T)) lz Y^f=l ^3 ,t{.Ry,D{T)) 


m=2 
oo 


dj' 


+ 


dx 


®p| 5 X^( 1 ^ 


(3.32) 


Furthermore, under Condition A2, there exists a positive constant C such that 


Y [R 


l Y,D(T) 



£(T||xi - x 2 1|) £(T||x 2 - xi| 


1 


D JD 


(3.33) 


C 



C(T) 

1 


D JD ll x l -X 2 


C(T) ||xi-x 2 | 
-^dxidx 2 = CTV (/C 2 ) < 00 


12 a 


dx 3 dx 2 


4-Tr ( R 


d; 


l Y,D(T) 


1 


£(T||xi - x 2 ||) £(T||x 2 - x 3 ||) C(T ||x m - xi| 


< C 


[C(T)] m J D (m) J D ||x 1 -X 2 || a 

1 


X 2 - X 3 | 


|X m - Xi 


dx 1 ... dx ?1 


Id ( m ) jb xi — 


X 2 || a ||x 2 - x 3 ||“ 


X m - Xi 


dx 3 ... dx r , 


(3.34) = CTr (/C™) < 00 , m> 2, 

since ||/C™||i < M||/C 2 ||i, for m > 2 and M > 0. 

From equations (I3.32D - (|3.34I) . for every T > 0, 

1 z rii (1 ( — l) m 1 /" 2z 

Wt( z )\ = exp - 


(3.35) 


(3.36) 


(3.37) 


exp 


exp 


cos 


+ sin 


exp 


exp 


exp 


2 ' 2 m — 2 \dx 

m =2 v 1 


2m—2 


Tr R 


1 (-l) m / 2^\ 2m_3 

2 ^ 2m-3 \d^ J 

171=0 


Y2m—2 
r ,D(T) 


Tr ( R 2m_3 
11 1 n Y,D(T) 


OO 

iX 


2m—2 


2 ^ 2m — 2 \d,T J 

m=2 v 7 

°o / -i\m / o.. \ 2m—3 

2 ^ 2m — 3 V dy 7 

v m=3 


TO ft 


1 OO 

tX 


(-l) m f2z\ 


2m—3 


2 2m. — 3 \dx J 

\ m =3 7 

-1 yv (-l) m /2z\ 

~2~ E 2m - 2 IdW 

, ra =2 x 7 


Tr ft 


>2m—3 
l Y,D(T) 


2m—2 


Tr ft 


->2m—2 

V,D(T) 


-1 y, (-l) 2n /2z\ 

IT ^ 4n - 2 V dr J 

7=1 7 


17 

OO 


4n—2 




-1 (- 1 ) 


E 

17=1 


2n+l 


4n 


2z\ 

dr) 


4 n 

1 Tr ( ft 


->4n 

l Y,D(T) 


( - ir ~‘ ^.rv (ft--) 


2m—3 \ 
Y,D(T) J 


i 1 / 2 
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(3.38) 



where we have applied, inside the argument of the exponential, the straightforward identities 
j 2 m —2 _ ^ 2 yi-i _ ^ | eX p^ u )| = C os 2 (u) + sin 2 (u) = 1, and the fact that the sequence 

of natural numbers m = 2,3,4, 5, 6 • • • = N —{0,1} can be obtained as the union of the sequences 
{2m — 2} m >2 = 2,4,6,... and {2m — 3} m >3 = 3,5, 7,... . Hence, in the above equation, the sum 
in N —{0,1} can be splitted into the sums X^m= 2 (~l) m-1 /(2m — 2) and ^“ =3 (—l) m /(2m — 3). 
Moreover, in (13. 3711 . we consider the sequence {2m — 2} m >2 = 2,4, 6,8,10,12,... as the union 
of the sequences {4n — 2} ri >i = 2,6,10,... , and {4n}n>i = 4, 8,12,... , corresponding to the 
changes of variable m = 2 n and m = 2n + 1. Thus, for m = 2n, 2m — 2 = 4n — 2, and, for 
m = 2n + l, 2m —2 = 4n + 2 — 2 = 4n. The sum ^m= 2 ( — 2) can then be splitted into 
the two sums X^i( — l) 2n /(4n — 2) and Xm^=i(—l) 2n+1 /(4n). Furthermore, the last identity 
in (I3.38P is obtained from the Fredholm determinant formula 


(3.39) V K * = det (I - w/C£) = exp 


OO 


E 



= exp 


EDmO 


• CcT 


k= 11=1 


of /C 4 at point oj = 16z 4 , which is finite for |w| < ^ , i.e., for \z\ < 

Theorem |Tj (see also equations (13.11) and (13.231) 1. 


- .,, , since, from 

4||^||J /4 


(3.40) Tr (/C 2 ) = [ [ || 1 ipa ^y = < 00 ■ 

j d j t~) ii y ii 

Hence, ||/C 4 ||i < M||/C 2 ||i = MTr (/C 2 ) < oo, for certain M > 0. 

From (|3.38l) . there exists ip(z) = lim^-nx) \ij) T (z)\ < oo, for 0 < z < [l/16||/C 4 ||i] 1(/1 . An 
analytic continuation argument (see [26], Th. 7.1.1) guarantees that ^ defines the unique limit 
characteristic function for all real values of z. 

From (|3.38|> . we now prove that ip(z) = ^(z), with i/)(z) given in (13.30[) (13.311) . 

Consider, for each z€l, the sequence 


(3.41) 





m 

Y,D(T) 


^ m> 2 


involved in the series appearing in the argument of the complex exponential defining tp T (z), in 
equation (13.32p . Keeping in mind that ( ij 2m = (—l) m , and that (i) 2m+1 = (—1 ) m i, let us split 
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(13.41 [) into the subsequences: 


(3.42) 


(3.43) 


1 — 


(-l) m (2z\ 


2m \d,T J 


2m 


— TV R 


->2 m 

l Y,D(T) 


m> 1 


{Qm} 


m> 1 — 


2m+l 


( 1) * ( Xl- ( R 2m+1 


2m + 1 \d,T J 

Under Condition A2, the following two inequalities hold: 


m> 1 


(3.44) 


(3.45) 


2m 




2m \cIt J 


Y,D(T) 


< C^-(2\z\) 2m Tr(JCl™) 


2m+l 


( 1) * Xr fi? 2m+1 

11 1 n Y,D(T) 


< C- 


1 


2m + 1 


(2|z|) 2m+1 Tr (/C 2m+1 ) • 


2m + 1 \dx J 
Moreover, for \z\ < 1/2, 

OO _« OO /'"V 

(2|4) 2 -" T rK™) < ]T-( 2 | 2 | riV(C) 

([P k| (2| 2 |)]- C 


m=l 


(3.46) 


(3.47) 


ra=l 
= In 


< oo 


/I 00 

£ 5^ (2 w> 2m+1 44C” +1 ) < ^;-(2MrMiv(C) 


m=l 


m= 1 
= 111 


([V Kl (2\z\)\ 


-CM 


< OO, 


considering \z\ < | A where Definition |Tj of the Fredholm determinant of /C 2 has been 

applied, since, as commented before, /C 2 is in the trace class from Theorem [I] Here, M is a 
positive constant such that, for every m > 1, Tr (/C 2m+1 ) < MTr (/C 2m ) , in the case where 
\\Y\\c(lRd)) < 1- Otherwise, 

Tr (/C 2m+1 ) < M'Tr (/C 2m+2 ) 


< M' 


sup \k{ICa) 
k>i 


Tr (/C 2m ) . 


Hence M = M' [sup fc>1 Afc(/C 2 )] 2 . 

Inequalities (I3.46P and (13.4711 allow us to apply Dominated Convergence Theorem for inte¬ 
gration with respect to a counting measure obtaining, for \z\ < ^ A YFRjO) ’ 


lim > 

T —ion ■ ^ 


(-l) m /2z\ 


T^too ^ 2m \dr J 

m= l v 2 


2m 

I Tr ( /?; 


-y2m 


l Y,D(T) 


E lim 


(-l) m /22 


(3.48) 


(-l) m i /2z\ 2m+1 


lim . . 

T—>oo 2m + 1 \dT J 

m= 1 v 7 


Tr ( /? 2m+1 ^ 
ri i jy y , D {T) ) 


m= 1 
oo 

E 

m=l 

oo 


T->oo 2m ydr 
-1) 


2m 


Tr ( R 2 ymt)) 


2m 


(2 z) 2m Tr (/C 2m ) 


m= 

oo 


um EhE (E) 

' T—>oo 2m + 1 \dx J 


Y lim 

T- 


2m+l 


Tr ( R 2m+1 
11 1 n Y,D(T) 


E//vt (2z) 2m+1 tv (C ,+1 ) 


m=l 


(3.49) 




























NON-CENTRAL LIMIT THEOREMS 


15 


where we have considered that, from (11.31) . for every m > 2, 

Dm 

V,D(T) 


(3.50) 


JSo" d? 


= Tr (K %), 


which is finite from Theorem [1} that ensures the trace property of K? a . Hence, for every m > 2, 
and z, 

(3 - 51 > & m (£) Tr A A = ™ <*>" Tr = m <2izr 

with c m being given in equation (13.311) . 

From equations (13.321) . (|3.42j) . (|3.43 j) . ()3.48j) . ()3.49l) and ((3.511) . we obtain, for \z\ < ^ 2 TV(a: 2 


)’ 


/ 1 °o i /9'\ m 

= ^“pUE^I^t) Tl i R y.D(T) 


2 ' m \ dr J 

\ m =2 ' ' 

f 1 °° \ (l °° 

T^ exp 2 S exp 2 S ^ 

\ tv 1 / \ m=l 


= exp 


1 00 \ /l 00 

?E r p UEri ft 


(3.52) 


m= 1 

i A i 


m=l 


exp 9 “ (2H) m Tr (/C™) = 

\ Z m =2 ” / 


An analytic continuation argument (see [26], Th. 7.1.1) guarantees that ip defines the unique 
limit characteristic function for all real values of z. 


We now turn to the proof of (ii). Under Condition Al, since H(||x||) < 1, and -B(O) = 1, 
we have 

B 3 (||x||)<B 3 (||x||), j> 3. 

Hence, 

-i 2 


Kt 


i 


£2 (vj'yj'ld-la 


[ G(Y(x)) dx - CgT d \D\ - [ H 2 (Y{x)) dx 

JD{T) 1 JD(T) 


oo (rjH^2 


E 


ye(T)T^\ ji j D(T) j D(T) 


D(T) 

H J (||x- y||)dxdy < 


(3.53) 


< 


£ 2 (T)T M-2q 


Id(t) Jd(T) 


B 3 (||x-y||)<fx<i y 


“ (Cf f 


3 =3 


J- 


By Condition Al, for any e > 0, there exists Aq > 0, such that for||x — y || > A 0 ,H(||x-y||) < 
e. Let V 1 = {(x, y) € D(T) X D(T) : ||x - y|| < A 0 )}, Vo = {(x, y) € D(T) x D(T) : ||x - y|| > 

Ai)}, 


(3.54) f f H 3 (||x-y||)dxdy = { / [ + [ [ \ £ 3 (||x - y||)dxdy = sft + s!p 
Jd(t)Jd(t) U Jv i J Jv 2 ) 
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Using the bound -B 3 (||x — y||) < 1 on T>i, and the bound -B 3 (||x — y||) < eB 2 ( ||x — y||) on 
V 2 , we obtain, 


5, 


(i) 


< 


T> i 


|S 3 (||x-y||)|dxdy <M x T d 


for a suitable constant Mi > 0, and for T sufficiently large, under A2, 


5, 


( 2 ) 


< 


< e 


V 2 


|5 3 (||x-y||)|dxdy <e / / R 2 (||x - y||)dxdy 

I'D 2 


[ [ R 2 (||x-y||)dxdy = e 

JO(T) JD(T) 


= e£ 


(3.55) 
Thus, for 


—2 CX I I 

Jd( 1) Jl 

< eCC 2 (T)T 2d ~ 2a [ 

Jl 


D(T) JD(T) ll x “ y 

£ 2 Pl|x-y||) 




d(i)Jd(i) l|x - y\\ 2a C 2 (T) 
dxdy 


dxdy 


D( l) Jd( i) ll x ~ y 


1 2a 


< oo, 0 < a < d/2. 


M 2 = C 


f f 

Jd l) Jl 


dxdy 


D(i) Jd{ i) ll x — y 


12 a ’ 


from (I3.53IMI3.55I) . we have 

I\t C 

(3.56) < (Mi V M 2 ) 


1 

2 

^ (Cf ) 2 

f f 

C(T)T d ~ a 


1 

-u 
_1 

Jd(t) Jd(T) 


+ e 


S 3 (||x-y||)dxdy 

p'2d—2a£2^~p ^' 


£, 2 (T)T 2d ~ 2a C 2 (T)T2d—2ct 


is arbitrarily small together with e > 0 as T oo. The desired result on weak-convergence 
then follows. □ 


Remark 5. Note that the proof of Theorem\^fi) is slightly different from the proof of Theorem 
2 in f23], for the case of non-linear functionals of chi-squared random fields. Specifically, we use 
here similar arguments to those ones considered in that paper, to derive the existence of a limit. 
However, to obtain the explicit limit characteristic function, we apply Dominated Convergence 
Theorem in a different way, considering equations \3.fl\j-\3.f3\). 


Remark 6. Consider the case of d = 1 and discrete time. That is, let {Y(t), t € Z} be a 
stationary zero-mean Gaussian sequence with unit variance and covariance function of the form 


m = 


m 


\t\ a 


for 0 < a < 1/2. The proof of the weak convergence result in [35] and [I2j is based on the fol¬ 
lowing formula for the characteristic function of a quadratic form of strong-correlated Gaussian 
random variables: 

= exp{—i zTdf 1 } [det (It — Tizdf 1 Rt)] ^ 


(3.57) =exp(f;(2i zdfJ) k ^ 

[k=2 



E 


T—l 


ex p 


t =0 
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where 

T—l T -1 

(3.58) -lk S P R T = ^2 ' ■' ~ *2|)-B(|*2 - *3|) ■ • • B{\i k - ill), 

d T d T il=0 i k= 0 

with dr = T l ~ a £.(T), Rt = E[YY'], Y = (Y(0),... ,Y(T — 1))', SpRx denoting the trace of the 
matrix Rt , and It representing the identity matrix of size T (see p. 39 of [25]). One can get 
a direct extension of formulae flff.57| ) and \3.58\) to the stationary zero-mean Gaussian random 
process case in continuous time {Y(t), t € M} (see [24]). but for d > 2 direct extensions of 
\3.51\ ) and 113. 5 8\) are not available. The present paper addresses this problem by applying alter¬ 
native functional tools, like the Karhunen-Loeve expansion and Fredholm determinant formula, 
to overcome this difficulty of discretization of the multidimensional parameter space. Note that 
the Fredholm determinant formula appears in the definition of the characteristic functional of 
quadratic forms defined in terms of Hilbert-valued zero-mean Gaussian random variables (see, 
for example, Proposition 1.2.8 in m- 

Remark 7. Expanding around zero the characteristic function H3.30\) . we obtain the cumulants 
of random variable Sqo, that is, n\ = 0, and 

(3.59) K k = 2 k ~ 1 (k - 1 )!c fc , k > 2, 

where c k are defined as in equation 113.31)) . The derivation of explicit expressions for c k would 
lead to the computation of the moments or cumulants of the limit distribution. This aspect will 
constitute the subject of a subsequent paper. 


4. Infinite series representation and eigenvalues 


The representation of the Rosenblatt-type distribution as the sum of an infinite series of 
weighted independent chi-squared random variables is derived in this section. As in the classical 
case (see Proposition 2 in m ), this series expansion is obtained from the double Wiener-Ito 
stochastic integral representation of in the spectral domain (see Theorem [3]). Proposition |T] 
and Corollary [2] below establish the connection between the eigenvalues of operator lC a in m 
and the weights appearing in the series representation derived. 

The following condition will be required for the derivation of Theorem [3](ii) below. 


Condition A3. Suppose that Condition A1 holds, and there exists a spectral density /o(|| A||), 
A G R rf , being decreasing function for ||A|| € (0, e], with e > 0. 

If Condition A3 holds, from equation (11.41) . applying a Tauberian Theorem (see m, and 
Theorems 4 and 11 in [22]). 

(4.1) / 0 (||A||)~c(d,a)£(p 

rf d ~ a ) 

Here, c(d,a) = 2 a ^ d/ 2 T ^ is defined in <EI)- 

Condition A3 holds, in particular, for the correlation function (11.51) . with the isotropic 
spectral density 


I ot—d 


0 < a < d, 


(4.2) 


/o(l|A||) 



a i a 

22- 1 7T2 


+1 



Kd 

2 


sin ^7 arg ^1 + exp 


1 + u& exp ( 



d 

u 2 du, 


where K v (z) is the modified Bessel function of the second kind. By Corollary 3.10 in [25], the 
spectral density dMD satisfies m , with a = /?7 < d. 
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The zero-mean Gaussian random field Y with an absolutely continuous spectrum has the 
isonormal representation 

(4.3) Y(x.) = [ exp(i(A,x)) y/f 0 (\\X\\)Z(dX), 

JR d 

where Z is a complex white noise Gaussian random measure with Lebesgue control measure. 


Theorem 3. Let D be a regular bounded open domain. 
(i) For 0 < a < d/2, the following identities hold: 


(4.4) 


/ 


\K (Ai + A 2 , D)\" 


dX x dX 2 


\d—ot 


'auMpyi 2 

,V2\D\. 


h(a)] 2 Tr(JC 2 a ) 

\D\ 2 


< oo, 


where is defined in \3.23 1) . 7 (a) is introduced in equation \3.3\) . and I\ is the char¬ 
acteristic function of the uniform distribution over set D, given by 

(4.5) K (A, D) = J e^pD (x) <bc=±- J e i(A ’ x) dx = 

with associated probability density function pr> (x) = 1/ \D\ if x S D, and 0 otherwise. 
(ii) Assume that Conditions Al, A2, A3 hold. Then, the random variable S/o admits 
the following double Wiener-ltd stochastic integral representation: 

(4.6) Sco = |B|c(<f,a) f ' g(A,,A 2 ) Z (<W Z , 

|| Al || 2 HA 2 H 2 

where Z is a complex white noise Gaussian measure with Lebesgue control measure, and 
the notation f R2 d means that one does not integrate on the hyperdiagonals X\ = ±A 2 . 
Here, the kernel H is given by: 


(4.7) 


and c (d, a) 


H (Ai, A 2 ) = K (A x + A 2 , .D), 
r (^) _ 1 

7 r d / 2 2 °T(a/ 2 ) 7(«)' 


Remark 8 . Our goal in this paper is to focus of the case of Hermite rank m = 2, which has 
very special properties not shared by the higher orders, such as the existence of eigenvalues. We 
are aware of the extension to all Hermite ranks, as described, for example, in the more general 
and different approach presented in the monograph by Major (1981). 


Proof, (i) From equation (I3.8|) and the proof of Theorem [H 


IM « 2 a-d = 


1 


1 


Id j(d-2a) J D ||x-y 
a 2 , 1 


12 a 
00 


dydx 


2 q(d — 2 a) 7 (d — 2 a) 


SH \2(r 2 \ - Tr (^a) ^ 


j = 1 


7 (d — 2 a) 


(4.8) 

since is in the trace class. Therefore, 1 o belongs to the Hilbert space H. 2 a -d with the inner 
product introduced in equation (13.911 . From equation (13.81) . we then obtain 


2 7 (d - 2 a) 




\D\ 


DWH2 a -d 


(2ir) d J Rd 


\K(ui, D)\ 2 \\<jOi || _ ^ _2a ^di 


,u>i. 


It is well-known that the Fourier transform defines an automorphism on the Schwartz space, 
which, in particular, contains C(f{D). Thus, the Fourier transform of any function in the space 
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7 -L 2 a—d can be defined as the limit in the space 'H 2 a-d of the Fourier transforms of functions in 
Cfi°(D). Therefore, from equation (|3.7|1 with /(z) = \D\ 2 \K(z, D)\ 2 , 


27 (d — 2a) 


Hence, 

(4.9) 


\D\ 


\D \ 2 7(2a) f 

( 2 vr) d [ 7 (a )] 2 J R d 
\D\ 2 'y(2a ) f 


2 ||u;i||“ d+2a duJi 


\K[ui,D)Y 


( 2 vr) d [ 7 (a )] 2 J R 2 d 


\D\ 


\K(\ 1 + \ 2 ,D)\' 


.7(a) 


\K(X 1 + X 2 ,D)\- 


\<jJl — (JL>2 

dX±dX 2 

(I|Ai||||a 2 ||)^- 

dX\dX 2 


-d+a U/ ^ n -d+a d 


\U>2\ 


duj 1 


A-| II II An 


\d—a ’ 


since — = 1. Note that, we also have applied the fact that, from Remark [2], 

lfl*lfl(x)= [ lD(y)lL»(x + y)dy = / l D (x + y)dy <E L 2 (D) C H 2a -d, 
J m d J d 


since 


lr)(y)ln(x + y)dy 


2 

dx < |H h(D) (0)| 3 , 


where, as before, \Bb,(D)\ denotes the Lebesgue measure of the ball of center 0 and radius R{D ), 
with R(D ) being equal to two times the diameter of the regular bounded open set D containing 
the point 0. Hence, J-(1d *1d)(A) = \D\ 2 \K(X, D)\ 2 belongs to the space of Fourier transforms 
of functions in di 2 a-d■ Summarizing, equation (|4.9I) provides the finiteness of (14.41) . i.e., assertion 
(i) holds due to the trace property of /C 2 for the regular bounded domain D considered (see 
Theorem [T]). 


(ii) The proof of this part of Theorem 0 can be obtained as a particular case of Theorem 5 in 
[22] (see also Remark 6 in that paper). Note that convexity is not used in the proof of Theorem 
5 of [22]. An outline of the proof of Theorem 5 in [22] for the case of Hermite rank equal to two 
is now given. 

Under Conditions Al, A3 (see also 63])), 


(4.10) F(x) = ^l/ exp (i (x, A)) K (A, D(T)) f^ 2 (X)Z(dX), x € D(T). 

(2n) a J R d 

Using the self-similarity of Gaussian white noise, and the Ito formula (see, for example, H21; 
[28]). we obtain from equation (14.101) 


S T = = 


T d ~ a C{T) J D{T) 

c(d,a)\D(T)\ 
T d ~ a C(T) 

c (d, a) | D 


H 2 (Y (x))dx 


d T d ~ a C(T ) 


f K(X 1 + X 2 ,D(T)) I -L- n/ 0 1 / 2 (A,) ] Z(dX 1 )Z(dX 2 ) 

J R 2 d l C (d, a) AA j 

f" K{X 1 + X2,D)\-±-f[fl / \x j /T)\ Z (dX\) Z (dX 2 ) ■ 
J R 2 d l c (a, a) 1 


(4.11) 
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By the isometry property of multiple stochastic integrals 


(4.12) 
where 

(4.13) 


E 


St — c (d, a) \D\ f H{ Ai, 
J R 2d 


Z(d\ 1 )Z(d\ 2 ) 

d—Oi d—a 

II Ai II 2 || Ao II 2 


1 2 


/ 

J R 2 


|K (Ai + A 2 , D)\ 2 [c (d, a) \D\] 2 Qt(Xi, A 2 V 


dXidX 2 


A nd—a || \ Md— a ' 

ill ll A 2|| 


Qt(Ai, A 2 ) 


|^ i ||(d-a)/2 ||^ 2 ||(d-a)/2 

T d ~ a C(T)c (d, a) 


n /o 1/2 (A/r) 


i=i 



From equation (14.11) . under Condition A3, we obtain the pointwise convergence of Qt( Ai, A 2 ) 
to 0, as T —>• oo. By Lebesgue’s Dominated Convergence Theorem, the integral converges to 
zero if there is some integrable function which dominates integrands for all T. This fact can be 
proved as in pp. 21-22 of [22], applying previous assertion (i) derived in this theorem. □ 


Alternatively, in the proof of Theorem (3[ii), the class CC of slowly varying functions, in¬ 
troduced in Definition 9 in m, can also be considered. Note that an infinitely differentiable 
function £(•) belongs to the class CC if 

1. for any 6 > 0, there exists Ao(5) > 0 such that X~ S C(X) is decreasing and X S C(X) is 
increasing if A > Ao(5); 

2. Cj E SC, for all j > 0, where £o(A) := C, Cj + i(A) := A£((A), with SC being the class 
of functions that are slowly varying at infinity and bounded on each finite interval. 

In that case, the following lemma can be applied in the proof of Theorem [3][ii). 

Lemma 3. Let a E (0 ,d), S E C'°°(s ri _i(l)), and C € CC. Let (£(x), x E be a mean-square 
continuous homogeneous random field with zero mean. Let the field £ (x) has the spectral density 
/o(u), u E M d , which is infinitely differentiable for all u/0 .If the covariance function B(x), 
x E M d , of the field has the following behavior 

(a) ||x||“B( x ) ~ S’(u^ji) £(||x||), x —s>oo, 
the spectral density satisfies the condition 

(b) ||u|| d - a / 0 (u) ~ S a4 C (ijAjy) , ||u|| — > 0. 

On the other hand, from Theorems|T]and[3](i), the spectral asymptotics of IC a and the Dirichlet 
Laplacian operator on L 2 {D) can be applied to verifying the finiteness of (|4.4|) for a wide class of 
domains D. Drum and fractal drum are two families of well-known regular compact sets where 
Weyl’s classical theorem on the asymptotic behavior of the eigenvalues has been extended (see, 
for example, [16] : [20]: [46] 1. In particular, as illustration of Theorem [3][i), we now refer to the 
case of regular compact domains constructed from the finite union of convex compact sets like 
balls, or by their difference which is the case, for instance, of circular rings. 


Examples. Let D = £q(0) U£>i((2, 0)) C R 2 , with £>i(0) = {( x\,x 2 ) E M 2 : \Jx\ + x\ < 1}, and 
£>i((2,0)) = {(xi,x 2 ) E K 2 : \J(x\ — 2) 2 + x\ < 1}- It is well-known (see [17], p. 57, Lemma 
2.1.3) that, for £>i(0) C M 2 and 0 < a < 1, 


Tr([/C® l(0) ] 2 ) = 


1 


/t?i(o) Jt3i( o) ll x - y 


1 2a 


dydx 


2 2-2a+l 7r 2-I r( 2^ ± l ) 
(2 — 2a)T(2 — a + 1) ’ 


where, to avoid confusion, for a subset S, we have used the notation /Cf to represent operator 
IC a acting on the space L 2 (S), and [/Cf] 2 = /Cf/Cf. 
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Hence, 


/ / 

Jb 1 COluBi (Y2,0l) Jb 


1 


Bi(0)UBi((2,0)) J Bi(0)UBi((2,0)) ll x ~ Y 


12a 


dyd'x. 


< 


/ / 

Jb 3 (o) Jb 


1 


b 3 ( 0 ) Jb 3 ( O) ll x - y 


12a 


dydx 


= Tr 


= 3 4_2 "Tr 


34-2002-20+1^2-^^ 2-20+1 ^ 
(2 - 2a)T(2 - a + 1) 


< 00 . 


(4.14) 


From Theorem (3)(i), equation (14.141) provides the hniteness of (14.41) for non-convex compact set 
£ = Bi(0)UBi((2,0)). 

These computations can be easily extended to the hnite union of balls with the same or with 
different radius, and to the case d > 2, considering the value of the integral 


/ / 

Jb r ( 0 ) Jb, 


1 


12a 


dydx = R 


2d— 2a. 


a rf (Si(0))i 


ib r ( 0 ) Jb r { 0 ) ll x - yii 

where the constant ad(B i(0)) is defined in (I3.26p . for 0 < a < d/2 (see [17], p. 57, Lemma 
2.1.3). 

For the case of a circular ring, that is, for 

D = B Rl (0 )\Br 2 (0) = {x E R“ : R 2 < ||x|| < 7?i}, R\ > R 2 > 0, 
we can proceed in a similar way to the above-considered example. Specifically, 

1 


< 


/ / 

Jb Ri (o)\b b jo) Jb, 

/ / 

Jbr. (o) Jb 


Ib Ri (o)\b R2 (o) Jb Ri (o)\Br 2 ( 0 ) ll x - y 


12a 


dydx 


Bfl, (0) Jb r . (0) ll x -yll 


= Tr ([^ l(0) ] 2 ) = i^" 2 "Tr ([/C® (0) ] 2 ) 


i? 4-2a 2 2- 2 a+ 1 7 r 2 -i r ( 2z^±i 


(2 — 2a)T(2 — a + 1) 


< 00 . 


From Theorem E)(i), equation (14.4j) is hnite for D = B Rl (0)\Bf{ 2 (0). Similarly, these computa¬ 
tions can be extended to the hnite union of circular rings. 


Remark 9. Note that for a ball D = £>i(0) = B( 0) = {x € ||x|| < 1}, the func¬ 
tion $(A) in H4-5\) is of the form exp (i (x, A)) dx = (2T) d / 2 ^^|jj^P , for d > 2, where 

J v {x) is the Bessel function of the first kind and order u > —1/2. For a rectangle, D = J/[ = 
{ai < Xi < hi, i = 1 ,..., d} , with 0 E f], we have i?(A) = n/ = 1 ( ex P (i Aj&j) — exp (iAjOj)) /iAy , 
for d > 1. Moreover for d = 2, considering the non-convex set D = B\(0) U £q((2,0)) C M 2 , 

i?(A) = i?(Ai, A 2 ) = f exp (i (x, A)) dx = 27r ^|| A ^ (1 + exp (2iAi)), 

J Bi(0)UBi(( 2,0)) ll A ll 

and for D = B Ri (0)\Br 2 (0), #(A) = (2 nRM (HAII^) /||A|| - (2ttR 2 )J 1 (||A||i? 2 ) /||A||. 


The following corollary is an extension of Proposition 2 in [12J. 

Corollary 1. Assume that the conditions of Theorem^ hold. Then, the limit random variable 
Soo admits the following series representation: 

















22 


N.N. LEONENKO, M.D. RUIZ-MEDINA, AND M.S. TAQQU 


(4.15) 

where 


See = c(d, a)\D\ Y T n (^)(4 - !) = Y X n( S °°)( £ l ~ !). 

nSN* nSNf 


c(d, a) = 


r(77 


ir d / 2 2 a r(a/2) 7 (a) 


was already introduced in \3.3\) . e n are independent and identically distributed standard Gaussian 
random variables, and n n (TL), n € N*, is a sequence of non-negative real numbers, which are 
the eigenvalues of the self-adjoint Hilbert-Schmidt operator 


(4.16) 


U{h){ X 1 )= / H 1 (X 1 -X 2 )h(X 2 )G a (d\ 2 ):L 2 E (R d ,G a ) 


T z 


\G a ), 


with 

(4.17) 


G a (d X ) = 


X 


I d—a 


dx, 


and L 2 E (j&L d ,G a ) denoting the collection of linear combinations, with real-valued coefficients, of 
complex-valued and Hermitian functions, that are square integrable with respect to G a (dx). Note 
that L|,(M rf , G a ) is a real Hilbert space, endowed with the scalar product 


(V’dV’ 2> Gq = / ^l( x )^2 (x)G Q (dx) 


(see [33] . pp. 159-161, for the case of L 2 E (M., d/3) spaces). The symmetric kernel H\ (Ai — A 2 ) = 
H( Ai, A 2 ), is defined from H introduced in equation 0 , in terms of the characteristic function 
K given in equation ra- 

The proof can be seen in Appendix A. Indeed, it constitutes an extension of Proposition 2 in 


In the following proposition the explicit relationship between the eigenvalues of JC a and TL is 

_ ^2 

derived. Note that Theorem[3li) provides the equality between the traces of operators 

and TL 2 , with, as before, TL having kernel H(-, •) given in equation (14.71) . 


Proposition 1 . The operators A a '■ L 2 E (W 1 ,G 0 


L|(R rf , 


G a ) 


A*(/)( Ar) = c{d, a) / H x (A : - A 2 ) / (A 2 ) G a (dX 2 ), 


and \D\ 1 JC a : L 2 (D ) —> L 2 (D) have the same eigenvalues. Here, c(d,a ) was already intro¬ 
duced in \3.3\), H 1 (Ai — A 2 ) = H( Ai, A 2 ) with kernel H being given in equation {» G a is 
introduced in \f.l r P) , and IC a is defined in 13.1) . 

The proof of this result is also given in Appendix A. (See [3H], for the case d = 1). 


Corollary 2. Let (A^Soo), k > 1} be the eigenvalues appearing in representation \f,15\ ), 
arranged into a decreasing order of the magnitudes of their modulus. Then, Theorem Q] holds 
for this system of eigenvalues. 


The proof directly follows from Corollary [fl Proposition Q] and Theorem [T] 
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5. Properties of Rosenblatt-type distribution 

This section provides the Levy-Khintchine representation of the limit random variable S 00 
(see @8], for d = 1, in the discrete time case), as well as its membership to a subclass of 
selfdecomposable distributions, given by the Thorin class. The absolute continuity of the law 
of Soo, and the boundedness of its probability density is then obtained. 

It is well-known that the distribution of a random variable X is infinitely divisible if for any 
integer n > 1, there exist X^ n \ j = 1,2,... ,n, independent and identically distributed (i.i.d.) 

random variables such that X = x[ n> + • • • + X^' 1 . Let TT>(M) be the class of infinitely divisible 

d 

distributions or random variables. Recall that the cumulant function of an infinitely divisible 
random variable X admits the Levy-Khintchine representation 


b 

(5.1) log E [exp (\0X)] = \ad — -9 2 T / (exp(i 6 ht) — 1 — ir(u)9)n(du), 9 € R, 

^ J — OO 

where a € R, b > 0, and 

(5.2) r(u) = 


u |it| < 1 

R M > 1 ’ 


and where the Levy measure /r is a Radon measure on R\ {0} such that /x({0}) = 0 and 

/ min(n 2 ,1 )/j,(du) < oo. 


An infinitely divisible random variable X (or its law) is selfdecomposable if its characteristic 
function (f)(9) = E[i6X\, 9 € R, has the property that for every c € (0,1) there exists a 
characteristic function cj) c such that (f)(9) = (f>(c9)(j) c (9), 9 € M. It is known (see [37], p.95, 
Corollary 15.11) that an infinitely divisible law is selfdecomposable if its Levy measure has a 
density q satisfying 

. , h(u) m 

q(u) = | . , uGR, 

|n| 

with h(u) being increasing on (—oo,0) and decreasing on (0,oo). Let ST>(M) be the class of 
selfdecomposable distributions or random variables. If Y £ ST>(M) then (see [19] ) 

n oo roo 

(5.3) Y= / exp (-s)dZ(s) = / exp(— s\)dZ(sX), A > 0, 

d J o d Jo 

where {Z(t), t > 0} is a Levy process whose law is determined by that of Y. 

We next define the Thorin class on M (see BU; 0 ; [18]) as follows: We refer to 'yx as an 
elementary gamma random variable if x is nonrandom non-zero vector in M, and 7 is a gamma 
random variable on R + . Then, the Thorin class on R (or the class of extended generalized 
gamma convolutions), denoted by T(R), is defined as the smallest class of distributions that 
contains all elementary gamma distributions on R, and is closed under convolution and weak 
convergence. It is known that T(R) C ST>( R) C TT>( R), and inclusions are strict. Since any 
selfdecomposable distribution on R is absolutely continuous (see, for instance, Example 27.8 
in m) and is unimodal (by [50]; see also Theorem 53.1 in [37]), then, any selfdecomposable 
distribution has a bounded density function. 

If a probability distribution function F belongs to T(R), then, its characteristic function has 
the form (see m, 0 ) 



+ 


iu9 

l+u 2 


U(du) , 


(5.4) 


(f)(9) = exp ( i 9a 
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where a € M, b > 0, and U(du ) is a non-decreasing measure on 
such that 


, called Thorin measure, 


f -1 


U(0) = 0, 


1-1 


|log |it|| U(du) < oo, 


1 


U(du) + 




;U(dll) < OO. 


The Levy density of a distribution from Thorin class is such that 

/o°° exp (-yu)U(dy), u> 0 


(5.5) 


\u\q(u) = 


Jo 00 exp (yu)U ( dy ), u < 0, 


where U(du ) is the Thorin measure. In other words, the Levy density is of the form /i(|tt|)/|it|, 
where h(\u\) = ho(r), r > 0, is a completely monotone function over (0, oo). 

The following result establishes the Levy-Khintchine representation of Soo, as well as the 
asymptotic orders at zero and at infinity of its associated Levy density. The membership to the 
Thorin self-decomposable subclass is then obtained. As a direct consequence, we then have the 
existence and boundedness of the probability density of Soo (see, for instance, Example 27.8 in 
[37]). 

Theorem 4. Let Soo be given as in Theorem {E with 0 < a < d/2. Let us consider Xk(Soo), 
k > 1, the sequence of eigenvalues introduced in Corollary U\ satisfying the properties stated in 
TheoremUl (see Corollary^. Then, 

(i) Soo € ZD(R) with the following Levy-Khintchine representation: 

(5.6) 


(5.7) 


oo, 


(f(9) = E[i9Soo] = exp j (exp(i uO) — 1 — i uQ) fi Q / d (du] 

where /a a / d is supported on (0, oo) having density 

^ OO 

q a /d(u) = 


k =1 


2 A k (S 0 


Furthermore, q a / d has the following asymptotics as u - 

[c(d,a)|D|( d -")/ d ] 1 /( 1 - a / d )r 


q a /d(u) 

a/d 


u > 0. 

-A 0+ and u — 
^-VT-a/d) 


2u[(\ — a/d)\ 




(a/d)-2 


(5.8) 


q a /d(u ) 


2 u 


[(1 - a/d)} 

exp(— u/2Xi(Soc)), as u —> oo, 


as u 


0 + , 


where c(d, a) is defined as in equation \3.13\) . 

(ii) Soo £ ST>(R), and hence it has a bounded density. 

(iii) Soo € T(M), with Thorin measure given by 


U(dx) = 


i fxL 
2 ^ 2\ k (Soo) 
k =1 


where 5 a (x) is the Dirac delta-function at point a. 
(iv) Soo admits the integral representation 


(5.9) 


= 


POO / °° \ POO 

J exp(-u)d\Y^Xk(So 0 )A (k \u)J=j exp (~u)dZ(u), 
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where 

OO 

(5.10) Z(t) = Y, \k(Soo)A^ k \t), t > 0, 

k= 1 

with A^ k \ k > 1, being independent copies of a Levy process. 


Proof, (i) The proof follows from Theorem [T] equation (13.121) . Corollary [2] and Lemma [4] below 
(see Appendix B), in a similar way to Theorem 4.2 in [38]. Specifically, let us first consider a 
truncated version of the random series representation (14.151) 

M 

s^ = Y x ^ s oo)(4-i), 

k =1 


with S™ —Soo, as M tends to infinity. From the Levy-Khintchine representation of the 
d 

chi-square distribution (see, for instance, [5], Example 1.3.22), 


E 



(5.H) 


M 

n^[exp(WA fc (5oo)(4-l))] 

k= 1 


M 

n«p 

k =1 
M 


+ /"(e Xp (i«„) -1) ( exp 
Jo 

exp (j (exp(i Qu) — 1 — i 9u) 


k =l 
exp 


2 u 

exp(-u/2A fc (S' 00 )) 


poo 

/ (exp(i#w) — 1 — i Qu) 

Jo 


2 u 


^ G x$/d) (exp(-u/2)) 


du , 



where G^^Ax) = \ To apply the Dominated Convergence Theorem, the 

following upper bound is used: 


(exp(i6Li) 


1 — i Qu) 


2 u U H°:/d) 


(exp(—n/2)) 


(5.12) 


^ ( ex P( — ' u /2)) 

q2 

< (exp(—u/2)), 


2 

where, as indicated in ]38], we have applied the inequality | exp(iz) — 1 — z\ < 4*-, for z € R. The 
right-hand side of (15.121) is continuous, for 0 < u < oo, and from Theorem (T] equation (13.121) . 
Corollary [2] and Lemma |4] in Appendix B, we obtain 


(5.13) 


uG Ka/d) (exp(—u/2)) 
u G\( a /d) (exp(-u/2)) 


~ uexp(—u/2Ai(5oo)), as u — > oo 

r (lVTTd) (1 - eM-WVr 1/(1 - a/d) 

a/d 

~ Cu as u — > 0, 


for some constant C. Since 0 < Ya/d < 1> equation ([5.131) implies that the right-hand side 
of (15.121) . which does not depend on M, is integrable on (0, oo). Hence, by the Dominated 
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Convergence Theorem, as M —? oo, 


E 


exr 


.p(i6>S'^ f) )l —> E [exp^Soo)] 


(5.14) 


= exp 


(exp(i0u) — 1 — i 6u) 


^■Gx( a /d) (exp(-u/2)) 


du 


which proves that equations (|5.6|i and (|5.7|) hold. 

Again, from Theorem [I] equation (1 3.12 [) . Corollary [2J and Lemma 2] below, 


2 u 


G x{a/d) (exp(—u/2)) ~ [h(d,a)|Z?| 1 -“/' i ] 1/(1 - Q/d); 


1 


1—at/cL 


u \-l/{l-a/d) 


(?) 


2u[(l — a/d)\ 


ct / d iii / \ (°i/d) — 2 

2 Wd[d(d,a)|.D| 1 - a / d ]i /(1_a/d) r (jrrTTTrf) u(W3) 


as u 


[(1 - a/d)] 

(5T5) —G A(a/d) (exp(-u/2)) ~ — exp(-u/2Ai(5oo)) as u 

Thus, equation (15.151) provides the asymptotic orders given in (15.81) . 


oo. 


(ii) From (i), it follows that S/o E <SX>(M), and hence it has a bounded density (see [8], Example 
27.8 in m, an d [50]). Note that an alternative proof of the boundedness of the probability 
density of S/o is provided in Appendix C, where an upper bound is also obtained. 


(iii) In view of (|5.5|) and (15.71) . Sqo E T(M) with Thorin measure given by 


(5.16) 


U(dx) = 2^ S ^l {x) ’ 

k= 1 


where 5 a (x) is the Dirac delta-function at point a. From Theorem [TJ Corollaries Q] and Propo¬ 
sition |T] (see also Corollary [2]) , the number of terms in the sum (|5.16p is infinite. Hence, the 
Thorin measure U(dx), as a counting measure, has infinite total mass. The form of Thorin 
measure is a direct consequence of (15.51) and (|5.7[) . 


(iv) As in m, we consider a gamma subordinator t > 0, with parameter A > 0, that is, 

a Levy process such that y A (0) = 0, and -P{ 7 a(^) € dx} = A _i r _1 (f) exp(— x\)x t ~ 1 dx, x > 0, 
and a homogeneous Poisson process N(t), t > 0, with unit rate. Assume that the two processes 
are independent. Then (see 0), for any c > 0, and A > 0, the representation (15.31) (see m) 
can be specified as follows: 

roc 

7a(c) = / exp(-t)d 7 A (JV(ct)). 
d Jo 

The process A(t) = r y 1 / 2 (N(t/2)) — t, t > 0, is a Levy process. 

For k > 1, let us consider 7 ^ (i) and A^ k \t) to be independent copies of 71 (i) and A(t), 

2 2 

respectively. Then, we have 

4 ~ 1 = 7i fc) Q) = exp (-u) dA {k \u), 

where are independent and identically distributed standard normal random variables as given 
in the series expansion (14.151) . Then, for A k{Soo), k > 1, being the eigenvalues appearing in 
such a series expansion, arranged into a decreasing order of their magnitudes, we obtain that 
the distribution of S ^ admits the integral representation (15.91) . with, A k , k > 1, in equation 
(15.101) being independent copies of the Levy process A(t) = 7 i/ 2 ( 77 (f/ 2 )) — t, t > 0. 

□ 
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For any 0 < a/d < 1/2, the Levy measure P’a/d satisfies 

/»oo 

/ u 2 n a/d {du) = E[S£j = [a d (D)} 2 . 

Jo 

Furthermore, when a/d —> 1/2, since (expert) — 1 — i 6u) —> (—1/2 )9 2 (see [HJ), we have 

( exp(i Ou) - 1 - \9u 2 f , A 

<P\P) = ex P I y - ^2 - M Va/dldu) 1 —> exp 

which means that Soo —> N( 0,1). 

In addition, from Theorem [H it can be proved, in a similar way to Corollary 4.3 and 4.4 in 
[3%j , that, for 0 < a/d < 1/2, the probability density function of S/o is infinitely differentiable 
with all derivatives tending to 0 as |x| —> oo. Also, the following inequality holds 

-P[Soo < —x] < exp > x > 0. 

We also note that, for e > 0, 

,. -P[Soo > u + e] _ f e \ 

™ P[S 00 >u] ~ eXP { 2A 1 (S 00 )J- 

Remark 10. In view of the integral representation \5.9\) . one can construct an Ornstein- 
Uhlenbeck type process 

dS(t ) = — XS(t) + dL(Xt), t > 0, A > 0, 

driven by a Levy process L(t), t > 0, and with marginal Rosenblatt distribution S oo. The driving 
process L(t ) is referred to as the background Levy process, and it is introduced in 15.101) . 

6. Appendices 
Appendix A 

Proof of Corollary [TJ. From condition (|4.4j) . operator TL is a Hilbert-Schmidt operator from 
L E (M. d , G a ) into L 2 E {R d , G a ), which admits a spectral decomposition, in terms of a sequence of 
eigenvalues {R n (TL), n € N*}, and a complete orthonormal system of eigenvectors {<^ n , n € N*} 
of L 2 E (JH d ,G a ), as follows: 



(6-1) tfi(x-y) = if(x,y) = Y 

neN$f 

where convergence holds in the L 2 E {R d ,G a ) <S> L 2 E {R d ,G a ) sense, i.e., 

2 


( 6 . 2 ) 


Y Ain(W)<Pn®Vr 

neN* 


= 0. 


L 2 E (R d ,G a )®L 2 E (R d ,G a ) 


We can establish the following isometry between the Hilbert space Lg(M rf , G a )<g>L 2 E (R d , G a ), 
and the two-Wiener chaos of the isonormal process X on H = L^j(M. d , G a ), given by 

(0 3 ) V . u ^ T'-i n ^ , vtu\ _ I ur„\ ^ (dx) 


X : he L 2 E (R d ,G c 


X{h) = [ h(x)- 

jRd 


|x||(^ a )/2 

(see [33], Chapter 9), considering the following identification between orthonormal bases of both 
spaces: For a given orthonormal basis {v? n <8>¥hL n G N*} of L^(M rf , G a )®L^(R d , G a ), its image 
by such an isometry I 2 is defined as 


(6.4) 


L 


MVn®Vn)= ^n( x l)^n( x 2) 


Z (dx 1 ) Z (dx 2 ) 

| Xl ||( d - a )/2 || x 2 || (^ —a )/2 ’ 
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which also defines an orthonormal basis in the two-Wiener chaos of the isonormal process X in 

m- 


From equations (16. 1 ji — (16.2 ji . by the orthonormality of the eigenvector basis <p n ®ip n of L 2 E (R d , G a ) 

(H, ® Vk) L 2 E (R d ,G a )®L 2 E {R d ,G a ) = ( Y ® Wn, <Ac ® W ) 

' nSN * ' L 2 ,(R d ,G a )®L 2 E (R d ,G a ) 

(6.5) = VkGNf. 

Again, from equations (16.11) (16.21) . and (16.51) . considering the isometry Z 2 in (16.41) 

!" HU if ' z < dXi > z{i ^ 

J R ™ 1 11 2 ||xi||( J -“)/ 2 ||x 2 ||(‘ , -“)/ 2 

r -'-—1 Z(rfXl) Z(rfX2) 

J R ; 


neNsf 


^n(xi)^n(x2) 


I11 a )/^ 11X 2 11 


( 6 . 6 ) 


<p n ( 

KjRd 


Z (dx) 


x 


d-a)/2 


where H 2 denotes, as before, the second Hermite polynomial. Note that summation and inte¬ 
gration can be swapped, in view of the convergence of the series (16.11) in the space L 2 E (R d , G a )<gi 
L 2 E (R d ,G a ), and the referred isometry between L 2 E {R. d , G a ) < 8 > L E (W l , G a ) and the two-Wiener 
chaos of isonormal process X introduced in (16.31) (see also equations (I6.2j) - (I6.5J) ). 

Note that 

f" , ^ Z ( dx ) 

n€N *’ 

are independent zero-mean Gaussian random variables with variance f R2 d |v3 n (x)| 2 G a ((ix), due 
to the orthogonality of the functions <p n , n E N*, in the space L E (R d , G a ). From equations (14.61) 
and (16.61) . 

Soo = c(d, a)\D\ Y M^)(4 - !)• 

a — 

neN* 

Equation (|4.15l) is then obtained by setting A n (S' 0O ) = c(d, a)|£)|// n ('H). 

Proof of Proposition |1] Let us consider T and J 7-1 the Fourier and inverse Fourier transforms 
respectively defined on L 1 (M rf ) and L 2 (R d ). Consider an eigenpair (//, h) of the operator A a , we 
have that J Rd |fi(y)| 2 ||y||d-Q < 00 . Applying the inverse Fourier transform T to both sides of the 
identity 

[ih = A a h, 

we get 

nX~\h) = X~\A a h) = c(d, a)J r ~ 1 (Hi * H 2 ), 

where, as before, 

Hi(Xi — A 2 ) = H(X 1 , A 2 ), 

with kernel H being defined in equation (14.71) . and iL 2 (y) = ||y|| _d " , " a fi(y). In the computation 
of this inverse Fourier transform, we note that H\ E L 1 (R d ) D L 2 (M rf ). In order to apply the 
convolution theorem, we first perform the following decomposition: 

H z(y) = llyir d+Q ^(y)ie 1 (o)(y) + llyir d+ “^(y) 1 R<i\B 1 (o)(y) : = # 2 ~(y) + - H 2 _ (y)> 

where £>i(0) denotes, as before, the ball with center zero and radius one in W l . Since 


[ h 2 (y)\\y\\- d + a dy<oo, 
JR d 












NON-CENTRAL LIMIT THEOREMS 


29 


H 2 € L 1 (M' i ), and € L 2 (M d ). Applying the linearity of the convolution and Fourier trans¬ 
form, the convolution theorem for both L 1 and L 2 functions (see [35], and [3D]) leads to 




(6.7) 


c(d, a)F~ l (H 1 * tf 2 ) = c(d, a) [F~ l {H 1 * # 2 ~) + -F _1 (77i * i7 2 + )] 
c{d,a)\D\- l l D (F- l {Hz + H+)) = 


where we have considered equations (14.51) and ()4.7D . From (16.71) . we can see that the support 
of J 7-1 )/;.) is contained in D, for any eigenfunction h of A a ■ The convolution theorem for 
generalized functions (see [45]) can be applied again to 77 2 , since h has compact support. By 
(14.171) . G Q (dx) = g a (x)dx, with g a fx) = ||x|| -d+ ". Then, 

My)lly|r d+Q =F(F-\h)*F- 1 {g a )) (y). 

Therefore, in equation (16.71) . we obtain 

tiF~\h) = c{d,a)\Dl l l D F- 1 [A{F- l {h)*F-\g a ))] 

(6.8) = c(d,a)|D| -1 l£) (j r_1 (/i) *-7 r_1 (fl , a )) . 


The inverse Fourier transform T 1 of g a ( y) = ||y|| d+a is obtained from equation 1 15.41 ) m 
Lemma[2](see also [41J, p. 117): 


-T (g Q )( z) = 


1 


7 r d / 2 2 a r(a/2) 


c(d,a)||z||“ T(^) 

Applying (]6.8() and this last relation, we finally obtain that, for an eigenpair (g, h) of A a , 


(6.9) gT 1 (h)( z) = \D\ x 1d(z) [ ||z - y|| a F 1 (h)(y)dy, 

Jd 

since, as commented before, J 7_1 (h) is supported on D. Thus, if (n,h) is an eigenpair of A ai 
then (g,F~ l {h)) is and eigenpair for \D\~ 1 IC a on L 2 (D). The converse assertion also holds. 
Hence, there exists a one-to-one correspondence between eigenpairs of A a and \D\~ 1 JC a , which 
preserves the eigenvalues, as we wanted to prove. 


Appendix B 

The proof of Theorem [4] is based on Lemma 4.1 in [38], which is now formulated. 


Lemma 4. Define the function G c (x) = YlT=i x ° k ^ w 'dh c = { c n } being a positive strictly 
increasing sequence such that c n ~ fin a , as n —> 00 , for some 1/2 < a < 1, and constant 
fi > 0. Then, 


( 6 . 10 ) 


G c (x) 

G c (x) 






x Cl , as x —> 0 

^ r (s) (1 " x) 


1/ot 

1 


as x 


1. 


Appendix C 

An alternative proof of the boundedness of the probability density of Soo, based on the series 
representation given in Corollary [T] is derived, and an upper bound is also provided. 
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Proof of boundedness of the probability density of S^. From Corollary [2l there exist 
two indexes ko and ki such that Ak 0 (S' oo ) > Ak^S'oo)- Then, 


5oo= E ^k(5oo)(4-l) =Ak 0 (S 00 )(4o-l) + Ak 1 ( 1 5 00 )(4 1 -l) + ^ 

keN^ 

where 

V= Y A k(Soo) (4 ~ 1) • 

keN^,k^k 0 ,ki 

Thus, 

Soo = Akj (Soo)(P^ 0 + 4j — (Ak 0 (‘S'oo) + Ak x (^oo)) + rj , 
where P = Ak 0 (5' 00 )/A kl (5 00 ). 

The random variables and £ t 1 are independent. Since the density of is of the form 



— t—^x 1/2 e x/2 , x > 0, 

T(i)42 


and the density of /?4 0 is given by 



- x > 0, 

PT{\)V2 K 


noting that ft = jg°°j > 1, then the density of ? = Pe^ n + 4, satisfies 


( 6 . 11 ) 


/?(«) = 


io 


/ej ( U - x)fp £ 2 o (x)dx 
l [ (u — x)~ 1 ^ 2 e^x~ 1 ^ 2 e~ 2 ^dx = 

o)v/5 Jo 


C kl 

e-“/2 

2r 2 ( 


2T 2 ( 


“TR n u - x) -W e 1(4),-!/!* 

o)v/5 ^0 


< 


■ 2 

e“/ 2 e 


, w i l 


2T 2 (±)v^ ./o 


[ (u — x) l / 2 x 1 / 2 da 

Jo 


< e 2 p 


Mij) 

2T 2 (i)v^ ^ 2^5 


< 


1 


< 


1 


Alc n (Soo) 2 
Ak 1 (‘S'oo) 


As the convolution of a bounded density with other is bounded, we then obtain the desired 
result. 
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